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Abstract of DE4423047 

The object reconstruction method involves projecting X-ray attenuation coeffts. along a straight line 
through the source and object to a detector is defined by an integral equation for different positions of 
the source. The three-dimensional structure of the object is reconstructed by means of the Feldkamp 
(1984) cone-beam algorithm with the region within the radiation cone parametrised into cylindrical co- 
ordinates. Three-dimensional back-projection is performed by addn. of filtered and geometrically 
weighted projection data or beam attenuation coeffts. and cyclic substitution of the co-ordinates of the 
points of intersection. 
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Priifungsantrag gem. § 44 PatG ist gestellt 

@ Verfahren zur Rekonstruktion eines Objekts 

(§) Beschrieben wird ein Verfahren zur Rekonstruktion eines 
Objekts, bei dem das Objekt von Strahlen, die von einer 
Strahlungsquelle mit einer kegelformigen Strahlgeometrie 
ausgehen, durchstrahlt und die Schwachung der Strahlen 
fiir verschiedene Positionen der Strahlungsquelle erfaSt 
wird, und bei dem der 3dimensionaIe Aufbau des Objekts 
mittels des Feldkamp-Atgorithmus rekonstruiert wird. 
Das erfindungsgemaBe Verfahren zeichnet sich daduroh 
aus, daS der innerhalb des Strahikegels liegende Bereich in 
Zylinderkoordinatan parametrisiert wird. 
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Beschreibung 



Die Erfindung bezieht sich auf ein Verfahren zur Rekonstruktion eines Objekts, bei dem Objekt von Strahlen, 
die von einer Strahiungsquelle mit einer kegelfdrmigen Strahlgeometrie ausgehen, durchstrahlt und die SchwS- 
5 chung der Strahlen fur verschiedene Positionen der Strahiungsquelle erfaBt wird, und bei dem der 3-dimensiona- 
le Aufbau des Objekts mittels des Feldkamp-Algorithmus rekonstruiert wird. 

Die Computertomographie hat sich beispielsweise in der medizinischen Diagnostik oder in der zerstdrungs- 
freien Materialpriifung sowie in der Geophysik bewahrt. Nur exemplarisch wird auf die nachfolgende Literatur 
verwiesen, auf die im iibrigen zur ErlSuterung aller hier nicht naher beschriebenen Begriffe ausdrUcklich 
to verwiesen wird: 

FUr medizinische Anwendungen wird verwiesen auf: 
Rontgen-Computer tomographic (CT): Louis et al, 1990, 
Single Photon Emission Computer Tomography (SPECT): 
Mangiosetal, 1989, 
15 Positron Emission Tomography (PET): Townsend et al. 1 989, 
Computer Angiography (CA) : Saint-Felix etal, 1993, 
Nuclear Magnet Resonance Tomography (NMR): Louis et al., 1 990, 

Fur Anwendungen bei der Materialpriifung wird verwiesen auf: 
Computertomographie (CT) mit Rontgenstrahlung: Maisl, 1992, 
20 CT mit Neutronenstrahlung: Maier-Schuler, 1992, 
CT mit -Strahlung: Kanamori et al„ 1985, 
CT mit Ionenstrahlung: Bench et al, 1991, 

Nuclear Magnet Resonance Tomography (NMR): Dieckmann et aL, 1991, 
optische Tomographic: Renz, 1992, 

25 seismische Tomographic ; Worthington, 1 984. 

Bei der zweidimensionalen Computertomographie wird eine einzelne Ebene im Objekt mittels vieler aus 
verschiedenen Richtungen aufgenommenen Projektionen rekonstruiert Urn jedoch die raumliche Lage, Form 
und Ausdehnung beispielsweise von Tumoren in der medizinischen Diagnostik bzw. von Rissen, Poren, Fremd- 
partikeln oder Legeirungsbestandteilen in der zerst6rungsfreien MaterialprQfung quantitativ ermitteln zu kttn- 

30 nen, ist die Untersuchung des gesamten Objektvolumens notwendig. Die Volumencharakterisierung durch 
Aufeinanderstapeln vieler Tomogramme mit geringer Schnitthtthe ist jedoch sehr zeitintensiv und daher fOr den 
praktischen Einsatz nicht geeignet. Abhilfe schafft hier die direkte dreidimensionaie Rekonstruktion von 
2D-Projektionsdaten, die unter vielen verschiedenen Winkein aufgenommen werden. Prinzipiell kann die dreidi- 
mensionaie Computertomographie durch folgende Verfahrensweisen reaiisiert werden: 

35 

a) parallele Durchstrahlung des Objektes (Fig. 1) 

b) Durchstrahlung des Objektes mit einer kegeiformigen Strahlgeometrie durch Verwendung einer punkt- 
formigen Strahlquelle (Fig. 2) 

40 Die parallele Durchstrahlung eines Objektes kann beispielsweise durch die Verwendung von Synchroton- 
strahlung [Flannery et al, 1987; Kinney et al., 1990) reaiisiert werden, Aufgrund der mangelnden Verfugbarkeit 
und der sehr hohen Betriebs- und Unterhaltskosten von Synchrotonstrahlquellen ist dieses Verfahren fur den 
kommerziellen Einsatz in der 3D-Computertomographie nicht geeignet. 

Kommerziel! verfugbare Tomographiesysteme mit einer punktfSrmigen Strahlenquelle Fuhrmann, 1992; 

45 Maisl, 1992) ermoglichen die DurchfUhrung der 3D-Computertomographie im Laborrahmen. Fig. 3 zeigt den 
schematischen Aufbau einer 3D-CT-Anlage mit Kegelstrahlgeometrie. Die Entwicklung von Rekonstruktionsal- 
gorithmen fur die 3D-Computertomographie unter Verwendung einer kegelfdrmigen Strahlgeometrie ist zur 
Zeit Gegenstand der Forschung. 
Die seit Ende der I970er Jahre entwickelten Computertomographiesysteme nutzen den kegelfdrmigen Strahl 

50 einer RSntgenquelle durch Verwendung zweidimensionaler Detektoren aus [Kowalski, 1979a; Robb, 1980]. 
Damit ergab sich die Notwendigkeit der Entwicklung neuer Rekonstruktionsalgorithmen, wenn der ausgenutzte 
Offnungswinkel des Strahlenkegels so groB ist, daB eine Annahme paralleler Strahlung nicht mehr gerechtfertigt 
ist. 

Zunachst wurden verschiedene Rekonstruktionsalgorithmen vorgeschlagen, die auf bekannten Algorithmen 
55 fur die zweidimensionale Rekonstruktion von Facherstrahlprojektionen basieren. Einige von ihnen erfordern, 
daB sich die Rontgenquelle (aber eine gedachte KugeloberflSche bewegt. [Nalcioglu et al, 1978; Denton et al, 
1979; Imiya et al, 1984; Peyrin, 1985]. 

Andere Rekonstruktionsverf ahren verlangen die Quellbewegung entlang eines Kreises [Feidkamp et al, 1 984] 
oder zweter paralleler Kreise [Schlindwein, 1978; Kowalski, 1979b] Martz et al. [199 1,1992 J Fuhrmann [1992] und 
60 Axelsson et al. [ 1 993]. 

Theoretische Arbeiten von, Kirillov [Kiriliov, 19611 Tuy [Tuy, 1983] und Smith [Smith, 1985] zur 3D-Compu- 
tertomographie mit kegelformiger Strahlgeometrie haben ergeben, daB eine exakte Rekonstruktion eines 
dreidimensionalen Objektes bei einer Bewegung der Quelle auf einer Kreisbahn um das Objekt nicht moglich ist. 

Welche praktisch realisierbaren Bahnen eine exakte Rekonstruktion ergeben (bei minimaler Redundanz der 
65 Quellpositionen) wurde von [Kudo et al, 19$0] diskutiert mit dem Ergebnis, daB zwei Bahngeometrien, namlich 
zwei aufeinander senkrecht stehende Kreise spiralfOrmige Umkreisungen des Objektes zu den besonders 
geeigneten Abtastgeometrien gehoren. 

Dariiberhinaus wurden von Tuy [Tuy, 1983] und Smith [Smith, 1985] theoretische Arbeiten zur Ldsung des 
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Rekonstruktionsproblemes fiir beliebige Quellbahnkurven durchgefiihrt Das von Smith [Smith, 1985] vorge- 
schlagene Verfahren wurde von Kudo [Kudo et al, 1990] und Smith [Smith, 1991, 1992] auf Computersystemen 
implementiert. Der Aigorithmus arbeitet mit der Hilberttransformation der ersten Ableitung der Radontransfor- 
mation des zu untersuchenden Objektes. Grangeat [Grangeat, 1991] hat ein 3D-Rekonstruktionsverfahren 
entwickeit, das auf der ersten Ableitung der Radontransformation des zu untersuchenden Objektes basiert. Aus 
dem Bereich der zweidimensionaien Computertomographie weiB man, daB Rekonstruktionsverfahren, die im 
Ortsfrequenzraum arbeiten, unter Umstanden weniger Rechenschritte benStigen als solche, die im Ortsraum 
operieren [Stark etal, 1981]. 

Martz et al. [1991,1992], Fuhrmann [1992] und Axelsson et ai. [1993] haben Grangeat-Aigonthmus modihziert, 
indem sie die Rekonstruktion nicht im Ortsraum sondern im Ortsfrequenzraum durchfiihren. Aktueller Gegen- 
stand der Forschung ist die Optimierung und Implementierung der verschiedenen Versionen des Grangeat-Al- 
gorithmus auf Parallel- und Vektorrechnern sowie der Entwicklung der dazu bendtigten speziellen Hardware- 
Architektur. , , , , , , „ , . 

Beziiglich der Anzahl der Rechenoperationen und der Qualit&t des Rekonstruktionsergebnisses tst der Feld- 
kamp-Algorithmus [Feldkamp et al., 1984] die bisher beste numerische L6sung des 3D-Rekonstruktionsprobie- 
mes. Praktische Erfahrungen zeigen, daS Rekonstruktionen mit dem Feldkamp-Algorithmus eine hohe Qualitat 
aufweisen, soiange der Ofmungswinkei des Rontgenkegelstrahls etwa 30° nicht uberschreitet [Webb, 1987; 
Kudo et al, 1990]. Ist der Offnungswinkei grdBer, so zeigen sich zunehmend Artefakte in den von der Zentralebe- 
ne weiter entfernten Schichten. In der Praxis wird hauf ig mit Offnungswinkeln nicht groBer als 20° gearbeitet, so 
daB der Aigorithmus sowohl einer guten Rechen- bzw. Speichereffizienz bzw. Rekonstruktionsqualitat derzeit 
am besten geniigt. 

Deshaib wird der Feldkamp-Algorithmus derzeit am haufigsten zur 3D-Rekonstruktion von tomographischen 
Daten eingesetzt [Webb et al, 1987; Bateman et al, 1989; Ohyama et al, 1989; Hu, 1989; Morton et al, 1990; Yan 
et al, 1991; Lin et al, 1992; Fuhrmann, 1992]. Die 3D-Rekonstruktionsverfahren von Smith [1991], Grangeat 
[1991] und Feldkamp et al. 1984] basieren auf dem Prinzip der gefilterten RUckprojektion [Brooks et al, 1976], 
Andere 3D-Rekonstruktionsverfahren, die z. B. auf der Matrixinversion, Reihenentwicklungsmethoden oder 
tterativen Methoden wie ART, SIRT etc. basieren [Censor, 1983; Chapman et al, 1986; Altschuier et al, 1977; 
Miggins et ai, 1 988] sind f Ur die 3D-Computertomographie nicht relevant, da die Rechenzeiten viei zu lange sind. 

Samtliche bekannte Verfahren haben jedoch den Nachteil, daB die Rekonstruktionszeiten far die 3D-Compu- 
terlomographie lang und haufig zu lang fttr den Einsatz in der Medizin oder der Industrie sind. 

Eine Verkiirzung der Rekonstruktionszeit ist prinzipiell mSglich durch 

a) Verbesserung der Rechnerarchitektur (Rechenleistung und Speicherkapazttat) 

b) Optimierung der Rekonstruktionsalgorithmen. 

Durch die fortschreitende Entwicklung der Computer stehen zwar immer leistungsfahigere Rechner zur 
VerfOgung, so daB Rechenleistung und Speicherkapazitat heutzutage kein grundsatzliches Problem mehr dar- 
stellen. Die Implementierung schneiler 3D-RekonstruktionsaIgorithmen ist aber trotz der fortlaufend erhdhten 
Leistungen der Computer nicht trivial und zur Zeit aktuelles Forschungsthe ma. 

Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren zur Rekonstruktion eines Objekts, bei dem das 
Objekl von Strahlen, die von einer Strahlungsquelle mit einer kegelfdrmigen Strahlgeometrie ausgehen, durch- 
strahlt und die Schwachung der Strahlen fiir verschiedene Positionen der Strahlungsquelle erfaBt wird, derart 
weiterzubilden, daB sich eine deutliche Verringerung der Rekonstruktionszeit und des Speicherbedarf gegen- 
iiber bekannten Verfahren ergibt. t 

Eine erfindungsgemaBe Ldsung dieser Aufgabe ist im Patentanspruch 1 angegeben. Weiterbildungen der 
Erfindung sind Gegenstand der abhangigen Ansprtiche. 

Der Erfindung liegt die Erkenntnis zugrunde, daB es moglich ist, bei der Rekonstruktion des 3-dsmensionafen 
Aufbaus des Objekts mittels des Feldkamp-Algorithmus die Rechenzeit dadurch deutlich zu verringern, daB der 
innerhafb des Strahlkegels liegende Bereich in Zylinderkoordinaten parametrisiert wird, und daB die 3D-RUck- 
projektion durch eine Addition der gefilterten und geometrisch gewichteten Projektionsdaten bzw. Strahl- 
schwachungskoeffizienten und eine zyklische Vertauschung der Schnittpunktskoordinaten durchgefiihrt wird 

Durch die erfindungsgemaBen MaBnahmen wird die Zahl der Rechenoperationen gegeniiber dem Feldkamp- 
Algorithmus Feldkamp et al, 1984] am ca. den Faktor 10 verringert. Die Zahl der Rechenoperationen, die beim 
Zylinderalgorithmus ftir die Schnittpunktsberechnung und die Interpolation vom Zylinderkoordinatensystem ins 
kartesische Koordinatensystem bendtigt wird, spielt f Ur die Gesamtrechenzeit keine Rolle, da sie eine GrSBen- 
ordnung kleiner (proportional N2) ist als die Zahl der Rechenoperationen, die f Ur die RUckprojektion (proportio- 
nal N3) benStigt wird. Der Speicherbedarf reduziert sich bei dem erfindungsgemaBen Verfahren urn ca. 30% 
gegeniiber dem Feldkamp-Algorithmus. Dennoch ist das Rekonstruktionsergebnis verglichen mit dem Ergebnis 
beim Feldkamp-Algorithmus vergieichbar. . 

Dies ermoglicht zusammen mit einem schnellen Computer Rekonstruktionszeiten, die einen breiten Einsatz 
der 3D-Computertomographie sowohl im medizinischen Bereich (CT, SPECT, PET, CA NMR) als auch im 
Bereich der zerstdrungsfreien PrQfverfahren (CT NMR, optische Tomographic) und in der Geophysik (seismi- 
sche Tomographic) erlauben. 

GemaB Anspruch 3 wird die Rekonstruktion der Strahlschwachungskoeffizienten |* mtt der Formel durchge- 
fiihrt: 
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hierbei bedeuten: 

M Anzahl der Projektionen, 

N Anzahl der Schichten in vertikaler Richtung, 

W Anzahl der WinkeKD 

NR Anzahl der Radialetemente, 

P gef ilterter Projektionswert an der Stelle (YS, ZS) 

GemaB Anspruch 4 ist es bevorzugt, daS das Verh&ltnis zwischen der Anzahl W der Winkel und der Anzahl M 
der Projektionen ganzzahlig ist. 

Weiterhin ist es bevorzugt, wenn fiir die Anzahl der Radialelemente gilt: 

NR = (N — l)/2. 

Weiterhin geniigt es, wenn die Berechnung der Schnittpunkte lediglich fiir den Bereich 
0° < 4>< 180° 

0 < R < n»radstep mit! < n < NR 
durchgefiihrt wird, 

Um die "gewohnte" Darstellung zu erhalten, wird gemaB Anspruch 7 das Rekonstruktionsergebnis aus dem 
Zylinderkoordinatensystem in em kartesisches Koordinatensystem transformiert. 

Die Erfindung wird nachstehend ohne Beschrankung des allgemeinen Erfindungsgedankens anhand von 
AusfUhrungsbeispieien unter Bezugnahme auf die Zeichnung exemplarisch beschrieben, auf die im iibrigen 
bezuglich der Offenbarung aller im Text nicht nilher erlauterten erfindungsgemaBen Einzelheiten ausdrQcklich 
verwiesen wird. Es zeigen: 

Fig. 1 eine schematische Darsteliung der 3D-Computertomographie mit paralleler Strahiung, 

Fig. 2 eine schematische Darstellung der 3D-Computertomographie mit kegelformiger Strahlgeometrie, 

Fig. 3 schematisch den Aufbau einer 3D-Computertomographie-Anlage, 

Fig. 4 eine schematische Darstellung zur Erlauterung des Feldkamp-Algorithmus, 

Fig. 5 ein FluBdiagramm fiir den Feldkamp-Algorithmus, 

Fig. 6 die Rekonstruktion eines dreidimensionalen Objekts beim Feldkamp-Algorithmus, 
Fig. 7 die Bestimmung der Intensitatswerte durch Interpolation, 
Fig. 8 die 2D-Computertomographie, 

Fig. 9 die 3D-Computertomographie mit kegelformiger Strahlgeometrie, 

Fig. 10 eine Erlauterung zur erfindungsgemaB vorgenommenen zyklischen Vertauschung, 

Fig. 1 1 die zyklische Vertauschung bei Drehung des Zylinderkoordinatensystems, 

Fig. 12 die Schnittpunktsberechnung, 

Fig. 1 3 die Abbildung eines Kreises, 

Fig. HdieSchnittpunktsellipse, 

Fig. 15 die Abbildung des Kreises auf etnen Detektorschirm, 
Fig. 16 die Veranderung der Schnittpunktsellipse, 

Fig. 17 eine schematische Darstellung zur effizienten Schnittpunktsberechnung, 

Fig. 18 die Interpolation vom Zylinderkoordinatensystem in das kartesische Koordinatensystem, 

Fig. 19 ein FluBdiagramm fiir das erfindungsgem&Be Verfahren, 

Fig. 20 ein FluBdiagramm fiir die erf indungsgemaBe Berechnung der Schnittpunkte, 

Fig. 21 die Ruckprojektion bei dem erfindungsgemaBen Verfahren. 

Im folgenden soil zunachst der Feldkamp-Algorithmus fiir die 3D-R6ntgen-Computertomographie beschrie- 
ben werden, wie der in Feldkamp et al, 1984 oder Rizo et aL, 1991 dargestellt ist 

Die Fig. 3 und 4 zeigen, daB die Projektion des Rontgenschwachungskoeffizienten u{x,y,z) iSngs einer Gera- 
den durch die Punkte S (Quelle) und M (Objekt) auf den Punkt A (Detektor) bei einer bestimmten Quellposition 
4> ist wie folgt definiert: 



Am Detektor A wird die Intensitat 1 gemaB dem Rontgenschwachungsgesetz fiir monoenergetische Strah- 
iung 




(1) 



(2) 



gemessen. I 0 ist die ungeschwachte Intensitat. 

Die Feldkamp-Rekonstruktion kann dann wie folgt beschrieben werden: 



f(M)=|-'jHD.yf(S > A)].g-dY 



(3) 



mit 



Yf(S.A) 




(4) 



wobei 



M Punkt im zu untersuchenden Objekt 

HD Entfaitung der Projektionsdaten YffSA) von der 1/r-Verschmierung 

0 Ursprung des Koordinatensystems 

O* Abbildung des Punktes O auf den Detektorschirm 

A Bildpunkt auf dem Detektorschirm 

S Position der Strahlenquelle 

O Winkel der Quellposition relativ zur x-Achse 

Die wesentlichen Schritte des Feldkamp-Aigorithmus sind in Fig. 5 dargestellt Die voflstandige Rekonstruk- 
tion eines 3D-Objektes erfordert die Wiederhoiung der in dem FiuBdiagramm in Fig. 5 dargestellten Rechen- 
schritte fur viete verschiedene Winkel (0° bis 360" ). 

Der Feldkamp-Aigorithmus ist eine Erweiterung der fur die 2D-Computertomographie gangigen Methode 
der geftlterten Riickprojektion [Maisl, 1992] auf 3 Dimensionen. Das Rekonstruktionsergebnis ist mit Ausnahme 
der in Fig. 6 dargesteiiten zentratebene, fur das es exakt ist, approximativ. (Fig. 6). Mit zunehmendem Abstand 
von der zentralebene treten aufgrund der Strahldivergenz Artefakte im Rekonstruktionsergebnis auf. 

Der Hauptrechenaufwand beim Feldkamp-Aigorithmus liegt nun in der dreidimensionaten Riickprojektion 
aller Einzelprojektionswerte. Dabei werden jeweils alle Volumenelemente adressiert und die entsprechenden 
gefilterten projekttonswerte dort abgelegt. Fur die Rekonstruktion eines Volumens mit N3 Elementen mit M 
Projektionen ist daher mit einem Rechenaufwand von ca. 10-M-N 3 Multiptikationen zu rechnen. Der Speicher- 
bedarf des Feldkamp-Aigorithmus liegt in der GroBenordnung 4«N 3 Bytes. 

Die dreidimensionale Riickprojektion setzt sich im einzelnen wie folgt zusammen. 



Fiir jedes Volumenelement muB der Schnittpunkt (YS,ZS) des entsprechenden Strahles mit dem Detektor 
berechnet und gespetchert werden (Fig. 6). Bei N 3 Volumenelementen ergeben sich dadurch 2-N 3 Schnittpunkt- 
skoordinaten. Die Berechnung der Schnittpunkte muB fiir insgesamt M Projektionen wiederholt werden. Zur 
Rekonstruktion eines Volumens mit N 3 Elementen und M Projektionen mit N 2 Daten ist also die Berechnung 
von 2»M'N 3 Schnittpunktskoordinaten notwendig. 



Fig. 7 zeigt, daB die Bestimmung der Intensitat an jedem Schnittpunkt durch biiineare Interpolation der 
Intensitatswerte an den benachbarten Detektorelementen erfolgt Jeder Schnittpunkt erfordert die Berechnung 
und Speicherung von 4 Gewichtungsfaktoren. Fiir M Projektionen und N 3 Volumenelemente miissen also 
4*M-N 3 Gewichtungsfaktoren berechnet werden. 

Die Fig. 8 und 9 zeigen, daB bei der 2D {3D)-Computertomographie mit Facherstrahl (Kegeistrahl-)geometne 
die Flache (der Raum) innerhalb eines Kretses (Zylinders) mit dem Radius Rmax rekonstrutert wird. 

Fur den Radius R msx ergibt sich: 

Rmax - adist * sin(atan(ydet/2bdist)) (5) ' 
hierbei bedeuten 

adist Abstand Quelie-Rotationsachse 



Schnittpunktsberechnung 



Biiineare Interpolation 



bdist Abstand Quelle-Detektor 

ydet Lange des Detektors in Y-Richtung. 

Voraussetzung fur die 2D- (3D)-Computertomographie ist, daB sich das zu untersuchende Objekt vollstandig 
5 innerhalb des Kreises (Zyiinders) mit dem Radius Rmax, das heiSt innerhalb des Strahlengangs, befindet. 
Die Erfindung geht nun von folgendem Grundgedanken aus : 

Der Bereich innerhalb des Kreises bzw. Zyiinders in den Fig. 8 bzw. 9 wird nicht wie ublich in kartesischen 
Koordinaten parametrisiert (Fig. 6), sondern in Polarkoordinaten (r f <D) bzw. in Zylinderkoordinaten (r,<D,z). 
Die Fig. 1 0 und 1 1 zeigen die Wahl der Koordinatensysteme, so daB der ubliche Zusammenhang gilt: 



r =^x 2 +y 2 
„ <p =atan(^) 

15 x 

z =z 

20 Bei einer Drehung des Polar- bzw. Zylinderkoordinatensystems urn den Winkel <&in C wird das Polar- bzw. 
Zylinderkoordinatenraster auf sich selber abgebildet. 

Aufgrund der Rotationsinvarianz des Polar-bzw. Zylinderkoordinatenrasters muB deshalb die Schnittpunkts- 
berechnung bei dem erfindungsgemaBen Verfahren nur einmal durchgeftihrt werden. Damit entfallt die Schnitt- 
punktsberechnung innerhalb der Riickprojektionsschleife in Fig. 5, wie sie beim Stand der Technik erforderlich 

25 ISt. 

Zunachst soli die Rekonstruktion filr die Zentraiebene (z =0) erlautert werden: 

Fig. 1 1 zeigt, daB bei der Drehung des Polarkoordtnatensystems um den Winkei <Di nc geht der Punkt (r,4>) in 
den Punkt (r,<D + Oi nc ) Ubergeht, 
Die Schnittpunkte (ys(r(i),O0))>zs(r(0.*(j) werden berechnet und in einer Tabeile abgelegt. Es muB gelten: 

30 

(p(j + 1 ) = <p(j) + (pine mit (pine = 2jt/M 

wobei M die Anzahl der Projektionen ist. 
Bei der Drehung um den Winkei 0>in C muB die Schnittpunktstabelle zyklisch vertauscht werden. Damit 
35 reduziert sich die Riickprojektion auf eine Addition der gefilterten und geometrisch gewichteten Projektionsda- 
ten und eine zyklische Vertauschung der Schnittpunktskoordinaten. 
Das Rekonstruktionsergebnis am Punkt u(r(i),<I>(j),z(k)) ergibt dann 

M W W 

40 ^(r(i)^(j),z(k)) = JP(YS(i,j +l.^,k),ZS(i,j +l^>k)) 

hierbei bedeuten: 

45 

M Anzahl der Projektionen 
WAnzahl der Winkei <D 

wobei P(YS,ZS) den gefilterten Projektionswert an der Stelle (YS,ZS) auf dem Detektorschirm darstelit Die 
50 Rekonstruktion des Rdntgenschwachungskoeffizienten u(r,<I>,z) fur das gesamte Volumen wird folgendermaBen 
beschrieben: 



mit 

eo NR Anzahl der Radialelemente (hier: NR = (N— 1)/2) 
N Anzahl der Schichten in vertikaier Richtung 

Gleichung 1 1 ist folgendermaBen zu verstehen: 

Die Summation der gefilterten und geometrisch gewichteten Projektionswerte P liber NR Radialelemente, W 
65 Winkelelemente und N Schichten in Z- Richtung entspricht der Riickprojektion der gefilterten und geometrisch 
gewichteten Projektionsdaten in das zu rekonstruierende Objekt fur eine bestimmte Projektionsrichtung. 

Dreht man das Zylinderkoordinatenraster um den Winkei <Dj nc (Fig. 11), so wird das Zylinderkoordinatenra- 
ster aufgrund der Rotationssymmetrie auf sich selber abgebildet Da die Schnittpunktskoordinaten nur einmal 
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berechnet und gespeichert werden mussen, reprasentiert die Summe liber N Projektionen in Gleichung 11 
sowohl die Verarbeitung aller M Projektionen mit jeweils N 2 Daten und die zyklische Vertauschung der 
Schnittpunktskoordinaten (YS,ZS). Damit die zyklische Vertauschbarkeit der Schnittpunktskoordmaten ge- 
wahrleistet ist, muB das Verhaltnis zwischen W (Anzahl der Winkei im Polar-bzw. Zylinderkoordmatenraster) 
und M (Anzahl der Projektionen) ganzzahlig sein. 
AufdiesenOberlegungenbautdieErfindungauf: 

In Fig. 12 ist der in Fig. 10 gezeigte Zyiinder in vier bezuglich der Koordinatenachsen symmetnsch liegende 
Bereiche unterteilt Aus Symmetriegrtinden mussen iediglich die Schnittpunkte der R6ntgenstrah!en fur den 
Bereich i, namlich 

0 ^ m • zstep < H mit 1 ^ m < (N- l)/2; 

0° £ <X> < 180° ; 0 < R £ n*radstep, 1 < n < NR 

berechnet werden, wobei zstep den Abstand der Schichten in z-Richtung und radstep den radialen Diskretisie- 
rungsabstand beschreibt 

Die Schnittpunkte flir die Bereiche 2-4 ergeben sich durch Symmetneiiberiegungen, das heiBt durch Spiege- 
lung an der Y*- bzw. Z*-Achse des Detektorkoordinatensysterns. Die Schnittpunktsberechnung ladt sich anhand 
der Abbildung eines Kreises mit dem Radius R=n«radstep (l<n^NR) und der Hdhenlage Z=m"zstep 
(1 < M < (N~~ l)/2) auf den Detektorschirm erklaren. Hierzu wird auf Fig. 13 verwiesen. 

Bildetmannuneinen Kreis in der H3he z « m • zstep und dem Radius R - n • radste p ab, so e rgibt sich auf dem 
Detektorschirm eine in Fig. 14 dargesteilte Schnittpunktsellipse, 

Die Schnittpunktsellipse kann allgemein durch folgenden Ansatz beschrieben werden: 



(Z±Z*) 2 ,(Y±Y M *)* , 
a 2 b 1 

Alle Koordinatenangaben, die mit einem (*) gekennzeichnet sind, beziehen sich auf das Koordinatensystem in 
der Detektorebene, alle anderen beziehen sich auf das Koordinatensystem im Objekt 

Das Zentrum der Schnittpunktsellipse (Fig. 13) ist (X'm, Y* M , Z*m). Die Haibachsen werden durch a und b 
beschrieben. Da der zu untersuchende Zyiinder (Fig. 9,10,12) symmetrisch bezilglich der X-Achse liegt, gilt: 



Y*m = 0. 

Aufgrund der Kegelstrahlgeometrie wird nicht der Mittelpunkt M eines Kreises (Fig. 13; dieser Punkt liegt auf 
der Rotationsachse in der H6he Z = m • zstep) auf den Mittelpunkt M e ff«(n) der Schnittpunktsellipse abgebildet, 
sondern der Punkt M e f<n) - (X(n),Y(n),Z(n)) (Fig. 14 und 15). 

Die Projektion der Strecke zwischen den Punkten M ef( <n) und P 3 (n) bzw. P 4 (n) (Fig. 15) auf den Detektor- 
schirm ergibt die Halbachse b der Schnittpunktsellipse (Fig. 14). 

Der Index m fur die Schnitthohe Z = m-zstep im Zyiinder (Fig. 1 2) ist folgendermaBen defimert: 



(l-N)/2<m£(N-l)/2. 

Fur m < 0 wird der Bereich X < X(n) (Fig. 15) auf den unteren Teil der Schnittpunktsellipse Z < Z*m (Fig. 14), 
derBereichX > X(n)(Fig. 15) jedoch auf den oberen Teil der Schnittpunktsellipse Z > Z*w(Fig. 14) abgebildet. 
Aufgrund von Symmetrietiberlegungen (Fig. 12) erhaltmanfur m > 0 folgendes Bild: 
Der Bereich X < X(n) (Fig. 15) wird auf den oberen Teil der Schnittpunktsellipse Z > Z* M (Fig. 14), der 
Bereich X > X(n) (Fig. 15) wird dagegen auf den unteren Teil m der Schnittpunktsellipse Z < Z*m (Fig. 14) 
abgebildet Punkte auf der Trennlinie X = X(n) werden auf die Halbachse b der Schnittpunktsellipse abgebildet. 

WahSt man m = 0, so entspricht dies der Abbildung der Zentralebene (Fig. 9) auf die Y*-Achse des Detektor- 
schirms (Fig. 14) mit Z* » 0. Symmetrieiiberlegungen in Fig. 12 zeigen, daB Iediglich die Schnittpunkte fur den 
Bereich 1 bestimmt werden mussen. 
Im folgenden soil die Berechnung der Halbachse a der Schnittpunktsellipse erlautert werden: 
Die GrdBe der Halbachse a ergibt sich durch die Abbildung der Strecke zwischen den Punkten M e ff(n) und 
P*(n) bzw. P 2 (n) auf den Detektorschirm (Fig. 14, 15): 

a 2 = m 2 • zstep 2 ■ (g a — gm) 2 

Die geometrische VergrSBerung am Punkt P,(n) bzw. P 2 (n) (Fig. 15) wird durch ga, am Punkt M ef f(n) durch g m 
beschrieben. Die Hohe der Z-Schicht wird durch Z= m • zstep angegeben, wobei zstep den Diskretisierungsab- 
stand in Z-Richtung beschreibt 

Die Berechnung der Halbachse b der Schnittpunktsellipse lauft wie folgtab: 

Die GrSBe der Halbachse b ergibt sich — wie den Fig. 14 und 15 zu entnehmen ist — durch die Abbildung der 
Strecke zwischen den Punkten M e r<n) und P 3 (n) bzw. P^n) auf den Detektorschirm. 



b 2 =g» 2 -n 2 -radstep 2 •(! .. " 2 -™fetep 2 

hHisr* 



Fig. 15 ist zu entnehmen, daB der Radius des n-ten Kreises durch n * radstep beschrieben wird, wobei radstep 
den radialen Diskretisierungsabstand beschreibt. bdist ist der Abstand zwischen der strahlenquelle und der 
Detektorebene. 

Mit diesen Berechnungsschritten kann man durch Einsetzen der Gleichungen 13 und 14 in Gleichung 12 unter 
10 der BerUcksichtigung, daB gilt 

Z*M = g m (n) * mzstep (15) 

und 

Y*=g b .Y (16) 
so erhatt man 

20 

Z^OO -m-au ! p.(g.(a).(-4)' +<*.<„) -g. („».(-!)*. II ^ SV * 

T 8. ( «) ' « ' radstep -(bdist 1 -tr • radstep ) 

25 miti = 1,2 und j = 1,2. 

Wahitmanj = 2, so wird die Schnittpunktseliipse far Z* > 0(m =4) berechnet. Wahit man anschIieBendj = l, 
so wird die bezOglich der Y*-Achse symmetrisch liegende Schnittpunktseliipse fUr Z* < 0(m = --4) berechnet 
Dies ist auch Fig. 14 zu entnehmen. 

Die Haibachse b der Schnittpunktsellipse teilt diese in zwei symmetrisch zueinander liegende Aste ein. 
30 Mit dem Index i wird nun festgelegt, welcher Ast der Schnittpunktseliipse berechnet wird (Fig. 14). 

Fur j = 2 und i = 2 wird der durch die Punkte P3*(n), P 2 *(n) und P 4 *(n) gekennzeichnete Teil der Schnittpunkt- 
seliipse fur Z* > 0 (m -> 4) berechnet. 

FOr j = 2 und i = 1 wird der durch die Punkte P3*(n), P t *(n) und P 4 *(n) gekennzeichnete Teii der Schnittpunkt- 
seliipse fur Z* > 0 (m =■ 4) berechnet. 
35 Fur j - 1 und i = 1 wird der durch die Punkte P 3 *(n), P 2 *(n) und P 4 *(n) gekennzeichnete Teil der Schnittpunkt- 
sellipse fUr Z* < 0(m =. -4) berechnet. 

FOr j = 1 und i » 2 wird der durch die Punkte P 3 *(n), P, *(n) und P 4 » gekennzeichnete Teil der Schnittpunkt- 
seliipse fUrZ* <0(m4) berechnet 

Die geometrische VergrSBerung eines beliebigen Punktes (X,Y,Z) im Kreis bzw. Zylinder (Fig, 10, 11) wird 
40 durch gb beschrieben. 

Da Gleichung 17 linear von m abhangt, kann man den Index m in Gleichung 17 durch folgende Beziehung 
eliminieren: 

AZ»g t „ = Z*ij, m+1 , n -Z*ij,m.n (18) 

45 

Man erhatt dann 



AZ^(Y)^tcp.( ga (u H ^ + < 8 ,( a )- g .(n)).(-4y.j; ZZZ^gZIZZ ) 

(19) 



55 Aus Gleichung 19 und Symmetriebetrachtungen gemaB Fig. 12 ergibt sich, daB 

a) die Schnittpunktskoordinaten Y* lediglich far die Punkte innerhalb des Halbkreises (0° :£ O < 180°, 0 :£ 
R Rmax, Z = 0) be rechnet werden mGssen. 

60 Fig. 1 7 veranschaulicht dies Symmetrieuberlegungen (Spiegelung an der Z*-Achse des Detektorkoordinaten- 
systems). Weiterhin ist zu beachten, daB die geometrische VergrSBerung lediglich far die Schnittpunkte in alien 
anderen Z-Ebenen unverandert bleibt (Fig. 16). 

b) die Schnittpunktskoordinaten Z* lediglich fUr die Ebene Z=m * zstep(m= 1,0° <, <S> < 180°, 0 < R < 
65 Rmax) berechnet werden mUssen. Die Schnittpunktskoordinaten Z* fur den Bereich 180° < flt> < 360° 

ergeben sich durch Symmetrieaberlegungen (Spiegelung an der Z*-Achse des Detektorkoordinatensy- 
stems,Fig. 17). Die Schnittpunktskoordinaten der Punkte innerhalb der Kreise mit (Z # 0,0° < O < 360°, 0 
< R < n»radstep mit 1 < n < NR) ergeben sich durch fortgesetzte Anwendung von Gleichung 18 und 19 
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bei entsprechender Wahl der Indices i, j, ra und n. 

Zur Schnittpunktsberechnung kann dann wie foigt vorgegangen werden: 
Gegeben ist em Zylinder mit folgender Diskretisierung: 

Radialrichtung: (N— i)/2 Eleraente, (N ganzzahlig und ungerade) 
Winkelrichtung: W Elemente, wo bei W ganzzahlig und gerade ist 
Z-Richtung: N Elemente (N ganzzahlig und ungerade) 
M: Anzahl der ftir die 3D-Rekonstruktion verwendeten Projektionen. 

Zur Berechnung der Schnittpunkte fur den Bereich 1 (Fig. 12) mit den Gieichungen 18 und 19 sind folgende 
Rechenschritte notwendig: 

Berechnung der geometrischen Verstarkung g m ftir die Punkte M e f<(n) = (X(n),Y(n) = 0, Z(n) = 0) ftir n = 
(N-l)/2Kreise(Fig. 15): 

(n) - bdist (20) 



x(n)=adist ^^_ m 

wobei bdist der Abstand zwischen der Strahlenquelle und dem Detektor und adist der Abstand zwischen der 
Strahlenquefle und der Rotationsachse ist (Fig. 13). 

Hierbei ist zu berticksichtigen, daB die Punkte M eft (n) - (X(n), Y(n)=0, Z(n) = 0) nicht notwendigerweise auf 
dem Zylinderkoordinatenraster liegen. Deshalb ist stets eine separate Berechnung und Speicherung von g m (n) 
notwendig. 

Die Berechnung der geometrischen VergrftBerung g a (n) fur die Punkte Pi(n) mit 1 < n < (N— 1)/2 (Fig. 15) 
erfolgt dann mit Gleichung 22: 

(22) 

Die X(n)-Koordinaten der Punkte Pi(n) ergeben sich aus Gleichung 6 bis 8. 

AnschlieBend erfolgt die Berechnung und Speicherung der geometrischen VergroBerung gb(X) ftir die Punkte 
P(RO)mit(0° < <D< 180°, 0 < R < n-radstepmit \ < n <, NRZ = 0)gemaB Gleichung 23: 

Die X-K.oordinaten der Punkte P(R3>) ergeben sich wiederum aus den Gieichungen 6 bis 8. 
Die weitere Vorgehensweise ist dann wie folgt: 

Berechnung der Y*-Koordinaten der Schnittpunkte fur die Punkte P(R,S>) mit (0° <; * < 180°, 0 <, R < 
n*radstep mit 1 < n <. NR,Z =0). 

Die X-Koordinaten der Punkte P(R,<J>) ergeben sich wiederum aus den Gieichungen 6 bis 8. 

Berechnung der d Z*-Werte fUr den Bereich (0° < <D < 180°, 0 < R < n«radstep mit 1 < n < NR, Z - m • 
zstep, m = 1) gemaB Gleichung 19. , 

Durch Anwendung von Gleichung 18 erhalt man die Z*-ICoordinaten der Schnittpunkte ftir Punkte in 
beliebigenZ-Ebenenmit(Z = m- zstep und (1 < m < (N — l)/2)sowie(0° < <!>< 180° ,0 <: R <, n-radstepmit 
1 < n < NR). 

Die Berechnung der Schnittpunkte fur die BereicheZ = m»zstep mit ((1 < m < (N-l)/2, 180° < <£>< 360°, 0 
< R £ n-radstep mit 1 £ n < NR)und((l~N)/2 ^ m £ ~1),0° S <t> <Z 360°, 0 < R <; n-radstep mit 1 < n < 
NR) erfolgt durch Symmetriebetrachtungen (Spiegelung an der Y*- bzw. Z*-Achse des Detektorkoordinatensy- 
stems.Fig. 12 und 14). 

Das erfindungsgemafie Verfahren iiefertdann — wie anschaulich Fig. 10 zu entnehmen ist — das Rekonstruk- 
tionsergebnis im Zylinderkoordinatenforrrfat Gangige 3D-Visuaiisierungssysteme erfordern jedoch, daB das 
rekonstruierte Volumen im kartesischen Gitterformat vorliegt (Fig. 6). 

Dies erfordert eine Transformation des rekonstruierten Volumens vom Zylinderkoordinatenformat ins karte- 
sische Gitterformat Fig. 18 zeigt fur eine beliebige Z-Schicht sowohl das Zyiinderkoordinatengitter als auch das 
kartesische Gitter. Zur Transformation des rekonstruierten Volumens vom Zylinderkoordinatenformat ins 
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kartesische Gitterformat bietet die Literatur Standardverfahren an [Stark et al, 1981; Matej et a!, 1990]. 

Unabhangig von der angewendeten Interpolationsmethode kann die Interpolation vom Zylinderkoordinaten- 
systera ins kartesische Koordinatensystem durch foigende Vorgehensweise optimiert werden. 
1st die Anzahl der Winkel W im Zyiinderkoordinatenraster eine gerade Zahl, so ist das Zylinderkoordinatenra- 
5 ster symmetrisch beziiglich der Koordinatenachsen X und Y. Dies bedeutet, da8 in Fig. 18 ledigiich die Interpo- 
lationskoeffizienten fur die Punkte in den Bereichen 1 und 2 (0° < «& < 90°) berechnet werden mflssen. Die 
Interpolationskoeffizienten in den Bereichen 3 bis 8 ergeben sich durch Symmetriebetrachtungen, das heiBt 
durch Spiegelung an den Koordinatenachsen. 
Ist die haibe Anzahl der Winkel W im Zyiinderkoordinatenraster (d. h. W/2) eine gerade Zahl, so ist das 
lo Zyiinderkoordinatenraster symmetrisch beziiglich der Koordinatenachsen X und Y und der Winkelhalbieren- 
den. Dies bedeutet, daB in Fig. 18 ledigiich die Interpolationskoeffizienten fUr die Punkte im Bereich 1 (0° < <D ^ 
45°) berechnet werden mussen. Die Interpolationskoeffizienten in den Bereichen 3 bis 8 ergeben sich durch 
Symmetriebetrachtungen, das heiBt durch Spiegelung an den Koordinatenachsen bzw. Winkelhalbierenden. 
Vorstehend ist die Erfindung anhand eines Ausflihrungsbeispiels beschrieben worden, Innerhalb des allgemei- 
t5 nen Erfindungsgedankens sind selbstverstandlich die verschiedensten Modifikationen moglich. In jedem Falle 
erhalt man aber durch die Reduzierung der 3D-Rttckprojektion (Fig. 19, 21) auf eine Addition der gefilterten und 
geometrisch gewichteten Projektionsdaten und durch eine zyklische Vertauschung der Schnittpunktskoordina- 
ten eine wesentliche Einsparung an Rechenzeit. 
Dariiberhinaus ist bet einer Weiterbildung der Erfindung eine besonders effiziente Transformation des Re- 
20 konstruktionsergebnisses vom Zyiinderkoordinatensystem (Fig. 10) ins kartesische Koordinatensystem (Fig. 18 
und22)m6g!ich. 

Nachstehend werden die Fundsteilen fur die vorstehenden Zitate angegeben. 
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PatentansprQche 

1. Verfahren zur Rekonstruktion eines Objekts, bei dem das Objekt von Strahien, die von einer Strahlungs- 
qudle mit einer kegelfdrmigen Strahlgeometrie ausgehen, durchstrahlt und die Schwiichung der Strahien 
fiir verschiedene Positionen der Strahiungsquelle erfaBt wird, und bei dem der 3-dimensionale Aufbau des 
Objekts mittels des Feldkamp-Algorithmus rekonstruiert wird, dadurch gekennzeichnet, daB der innerhaib 
des Strahlkegels iiegende Bereich in Zylinderkoordinaten parametrisiert wird. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB die 3D-Riickprojektion durch eine Addition der 
gefilterten und geometrisch gewichteten Projektionsdaten bzw. Strahlschwachungskoeffizienten und eine 
zyklische Vertauschung der Schnittpunktskoordinaten durchgefuhrt wird. 

3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, daB die Rekonstruktion der Strahlschwa- 
chungskoeffizienten \i mit der Formel durchgefuhrt wird: 



hierbei bedeutem 

M Anzahl der Projektionen, 

N Anzahi der Schichten in vertikaler Richtung, 

W Anzahl der Winkei 0> 

NR Anzahl der Radialelemente, 

P gefiiterter Projektionswert an der Stelle (YS, ZS). 

4. Verfahren nach Anspruch 3, dadurch gekennzeichnet, daB das Verhaltnis zwischen der Anzahl W der 
Winkei und der Anzahl M der Projektionen ganzzahlig ist 

5. Verfahren nach einem der Anspriiche 2 bis 4, dadurch gekennzeichnet, daB fiir die Anzahl der Radialele- 
mente gilt: 

NR =■ (N — l)/2 

6. Verfahren nach einem der Anspriiche 1 bis 5, dadurch gekennzeichnet, daB die Berechnung der Schnitt- 
punkte lediglich fur den Bereich 

0° < <D < 180° 

0 < R < n»radstep mit 1 < n ^ NR 
durchgefuhrt wird. 

7. Verfahren nach einem der Anspriiche 1 bis 6, dadurch gekennzeichnet, daB das Rekonstruktionsergebnis 
aus dem Zylinderkoordinatensystera in ein kartesisches Koordinatensystem transformiert wird. 



Hierzu21 Seite(n) Zeichnungen 
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